Shear Banding and Spatiotemporal Oscillations in Vortex Matter in Nanostructured 
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We propose a simple nanostructured pinning array geometry where a rich variety of complex 
vortex shear banding phenomena can be realized. A single row of pinning sites is removed from 
a square pinning array. Shear banding effects arise when vortex motion in the pin-free channel 
nucleates motion of vortices in the surrounding pinned regions, creating discrete steps in the vortex 
velocity profile away from the channel. Near the global depinning transition, the width of the band 
of moving vortices undergoes oscillations or fluctuations that can span the entire system. We use 
simulations to show that these effects should be observable in the transport properties of the system. 
Similar large oscillations and shear banding effects are known to occur for sheared complex fluids 
in which different dynamical phases coexist. 

PACS numbers: 74.25.Uv,74.25.Wx 
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Sheared complex fluids [H, S S 0, S Hi and colloidal as- 
semblies exhibit the dynamical phenomenon of shear 
banding, in which the motion of the system under an ap- 
plied shear becomes stratified or concentrated in certain 
regions. For many complex fluids where different dy- 
namical flow phases can coexist, coupling of these phases 
under shear can induce additional dynamics such as large 
oscillations and strongly fluctuating chaotic phases 0, Q ■ 
It would be interesting to explore whether similar behav- 
iors can occur in other systems such as magnetic vortices 
in type-II superconductors. Shear banding induced by 
the presence of periodic structures, such as a periodic 
pinning array for superconducting vortices, has not been 
considered previously. Shearing in complex fluids or col- 
loids in the presence of quenched disorder could be real- 
ized by having the flow traverse an optical trap array or 
a periodic array of micron-scale obstacles. 

In order to shear a superconducting vortex system, un- 
conventional contact geometries such as the Corbino ge- 
ometry are normally used (1, 0, [l^, 11 , l2| ■ Simulations 
of pin-free systems using such geometries show transi- 
tions from rigid body rotation to liquid-like flow states 
[iSs], and experiments revealed a similar effect near the 
vortex liquid to vortex solid transition In other ex- 
periments, noise and transport measurements of moving 
vortices were performed in a Corbino geometry IT^ . 
Recent simulations of Corbino contacts applied to very 
small superconducting disks produced a rich variety of 
novel vortex dynamics ^14 ] . In this work we propose the 
use of a simple periodic pinning array geometry rather 
than a Corbino geometry to shear the vortices. We ob- 
serve a rich variety of shear banding effects including an 
ordered phase with velocity steps, giant spatiotemporal 
oscillations, and strongly fluctuating chaotic phases. 

Using nanotechnology it is now possible to create tai- 
lored pinning array structures in which the lattice sym- 
metry and the order and size of the pinning sites can 
be carefully controlled 3, It has also been demon- 



strated experimentally that these arrays can be diluted 
by systematically removing individual pinning sites 17 1. 
Here we propose removing entire rows from a square pin- 
ning array and applying a magnetic field such that the 
vortex density matches the density of the original undi- 
luted square pinning array. Experimentally, this geome- 
try could also be achieved on a larger scale by periodically 
removing rows of pinning sites. We consider samples in 
which the pinning sites are sufficiently small that each 
pin can capture only a single vortex; thus, when there 
are more vortices than pinning sites, a portion of the 
vortices are located outside the pinning sites in the pin- 
free regions. Although this geometry is relatively simple, 
to our knowledge it has not been considered in previous 
simulations or experiments. Simulations of vortex dy- 
namics in square pinning arrays near and just above the 
first matching field showed that the interstitial vortices 
are more mobile than the pinned vortices and depin at 
lower driving currents, while coupling between the mov- 
ing interstitial vortices and the pinned vortices can pro- 
duce a series of transitions between different dynamical 
phases including a disordered or turbulent phase with 
a crossover to a laminar phase at higher drives [l^l- Al- 
though the simulations provided clear predictions for fea- 
tures in the transport curves associated with these tran- 
sitions, the experimental observation of these phases was 
not achieved until very recently due to new experimen- 
tal capabilities for reducing local heating effects and 
for performin g tr ansport measurements at temperatures 
well below Tc [20| . Molecular dynamics-style simulations 
should serve as excellent models for the experiments that 
are currently possible. 

In this work we demonstrate that when a row or rows 
of pinning sites are removed from a square pinning array, 
the vortex motion under an applied driving current first 
occurs in the pin-free channel and nucleates the depin- 
ning and motion of the vortices in the pinning sites. This 
produces a series of steps in the velocity response along 
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with fluctuating or chaotic phases and giant temporal 
oscillations when a band of moving vortices propagates 
away from and toward the pin-free channel. 

We simulate a two-dimensional system of supercon- 
ducting vortices interacting with a modified square pin- 
ning array. The sample is of size L x L and has periodic 
boundary conditions in the x and y directions. The undi- 
luted square pinning array contains Np pinning sites at 
a density of Up — Np/L'^, and the sample contains 
vortices at a density of n„ — N^/L"^ . The field at which 
Nv = Np is the matching field i?^. The dynamics of a 
vortex i at position is given by the overdamped equa- 
tion of motion 

n^ = Fr + + (1) 

The repulsive vortex-vortex interaction force is F™ = 
^^"j- Foifi(i?y /A)Rij. Here Ki is a modified Bessel 
function that decays exponentially at large distances, 
^0 = 0o/27r/ioA'^, A is the London penetration depth, 
4>o = h/2e is the flux quantum, /xg is the permeability of 
free space, Rij = jR^ — Rj|, and R^j ~ (R^ — Rj)/i?y . 
The pinning sites are modeled as parabolic traps with 
a maximum force of Fp and radius Rp — 0.35A which 
are placed in a square array with lattice constant a: 

= Ek' FpFoiR^k/RpMRp - R^k)'R^k, where 6 is 
the Heaviside step function. We set Up = 0.28/A^. We 
keep Rp and Fp small enough that a maximum of one 
vortex can be captured by a given pinning site, and we 
remove one row of pins. The initial vortex positions 
are obtained through simulated annealing, which pro- 
duces a square vortex arrangement at the matching fleld 
B — Bff,. The external driving force F'^^* = i^£)X rep- 
resents the Lorentz force from an applied current. We 
increase Fjj in small increments and measure the vortex 
velocity averaged over t — 8000 simulation time steps, 
(K) = N-W-^ dB.,{t)/dt • X. 

In Fig. 1(a) we plot the velocity (14) versus force 
Fo/Fp curve for a sample with Fp — 0.1 and B/B^ = 
1.0. The central portion of the sample is illustrated in 
Fig. 1(b) at Fo/Fp = 0.6, where only the interstitial vor- 
tices in the pin-free row are in motion. There are four 
distinct regimes in Fig. 1(a). In the pinned regime at 
low the interstitial vortices in the pin- free channel 
remain pinned due to the repulsive interactions from the 
vortices at the pinning sites. At higher drives, the in- 
terstitial vortices depin and enter the single channel flow 
(SCF) state illustrated in Fig. 1(b). A sharp increase 
in (14) occurs near F^/Fp = 0.75 and corresponds to 
the onset of motion of the vortices in pinning rows ad- 
jacent to the pin-free channel, termed the shear banding 
(SB) regime. The velocity of the vortices moving along 
the pinning rows is less than the velocity of the vortices 
in the pin-free channel, and during the time interval re- 
quired for the vortices in the pinning rows to move a 
distance a and hop from one pinning site to another. 
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FIG. 1: (a) The average velocity (Vx) vs applied force Fd / Fp 
for a system with Fp = 0.1, B/B^ = 1.0, and a single removed 
row of pinning. Four regimes are marked: a pinned regime, 
a single channel flow (SCF) regime, a shear banding (SB) 
regime, and a moving (M) phase in which all of the vortices 
are flowing. Inset: A blowup of the main panel for 0.081 < 
Fo/Fp < 0.0865 in the shear banding phase where a short 
time average is used for the velocities so that the fluctuations 
can be seen more clearly. Large scale oscillations occur as the 
system switches back and forth between three moving rows 
of vortices and all rows of vortices moving. On a long time 
average, the oscillating regime has an intermediate value of 
(Vx) as shown in the main panel, (b) The vortex positions 
(dots), trajectories (lines), and pinning site locations (open 
circles) for a small subsection of the sample in (a) in the SCF 
phase at Fo/Fp = 0.6. (c) Vx vs time for the system in (a) 
at Fo/Fp = 0.7, 0.85, and 0.92 (from bottom to top). 



the vortices in the pin-free channel move a distance 2a 
or greater. A small step in (14) appears when the mo- 
tion of the vortices in the adjacent pinning sites locks 
into step with the motion of the vortices in the pin-free 
channel. For Fo/Fp > 0.855 all the vortices are mobile 
and {Vx) increases linearly with increasing Fjj/Fp. For 
0.83 < Fn/Fp < 0.855, (14) takes a value that is close 
to half of the Ohmic value, indicating that on average 
half of the vortices in the sample are moving. It is not 
clear from the time-averaged data whether this means 
that half of the vortices are always moving and the other 
half are always pinned, or whether the number of mov- 
ing vortices is fluctuating. To clarify this, in the inset of 
Fig. 1(a) we plot (14) over a small range of Fo/Fp using 
a shorter time average of t = 175 simulation time steps. 
For 0.83 < Fo/Fp < 0.85, this (14) response exhibits 
large oscillations with a maximum velocity close to the 
Ohmic value and a minimum velocity close to that associ- 
ated with vortices flowing in three rows, indicating that 
the number of moving vortices is oscillating over time. 
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In Fig. 1(c) we plot the instantaneous Vx versus time at 
fixed Fo/Fp = 0.7, 0.85, and 0.92. For Fo/Fp = 0.7, 
only three vortex rows are moving and there is a high 
frequency periodic component of produced by the mo- 
tion of the vortices over the periodic pinning sites. For 
Fn/Fp = 0.92, all of the vortices are moving and there is 
again a high frequency periodic signal in Vx caused by the 
periodic motion over the pinning sites. At F^/Fp = 0.85, 
pronounced oscillations appear in Vx that have a much 
larger amplitude and a lower frequency, indicating that 
these oscillations are not associated with the periodicity 
of the pinning sites. A high frequency oscillatory com- 
ponent of Vx which is a result of the motion over the 
periodic pinning sites does appear; however, the larger 
scale oscillations are due to the emission and absorption 
of a transverse front of moving vortices originating at the 
pin-free channel. 

In Fig. 2 we illustrate the vortex trajectories dur- 
ing consecutive time intervals for the system shown in 
Fig. 1(c) at Fu/Fp = 0.85 for a single cycle of the large 
oscillations in Vx, which has a period of T = 15750 sim- 
ulation time steps. At the start of the cycle, Vx is at its 
minimum value and the vortices in the pin-free channel 
are moving along with vortices in the two adjacent pin- 
ning rows. Figure 2(a-c) shows that during the first half 
of the oscillation, t/T < 0.5, an increasing number of 
vortices in the pinning rows become mobile, and there is 
a well defined front between the moving and non-moving 
portions of the sample. At the halfway point of the oscil- 
lation, t/T = 0.5, Fig. 2(d) indicates that all of the vor- 
tices are mobile. When the two depinning fronts meet, an 
instability is triggered that causes the vortices to become 
pinned again, and two pinning fronts begin to propagate 
back toward the pin-free channel during the second half 
of the period, 0.5 < t/T < 1. Finally at t/T = 1, the 
system returns to the initial state containing only three 
moving rows of vortices. The period T of the oscilla- 
tions at fixed Fjj depends on the system size since the 
depinning and pinning fronts propagate at a fixed speed. 
By simulating large systems, we find that T scales with 
L. The high-frequency periodic component of the vortex 
motion that is generated when the vortices move over 
the periodic pinning array does not change with system 
size and is a function only of the pinning lattice constant 
a. Shear banding effects resulting in large oscillations 
have been observed for shearing in many complex fluids 
and are due to a nonlinear coupling between two differ- 
ent types of coexisting structures that form in the fluid 
[1, Q . In our vortex system the two structures are the 
square pinned vortex lattice and the distorted hexagonal 
ordering of the moving vortices. The mismatch between 
these two phases results in the formation of topological 
defects at the boundary between the phases. When the 
depinning fronts meet, it is possible that a portion of the 
defects annihilate, changing the direction of propagation 
of the front. 
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FIG. 2: The vortex positions (dots), trajectories (lines), and 
pinning site locations (circles) for the system in Fig. 1(c) at 
Fo/Fp — 0.85 during consecutive time intervals in a single 
large oscillation of Vx with period T. (a) At a minimum in Vx , 
t/T — 0, the vortices in the pin-free channel move along with 
vortices in one pinned row on each side of the channel, (b) 
At t/T — 0.125, a front of moving rows begins to propagate 
outwards from the pin-free channel and here three pinned rows 
on each side of the channel are moving, (c) Near t/T — 0.25 
the moving region continues to grow, (d) At t/T = 0.5 the 
front has moved completely through the sample and all of the 
vortices are moving, (e) Near t/T = 0.8 the moving region 
is contracting, (f) At t/T = 1 the system returns to the 
minimum state of three moving rows of vortices. 



In Fig. 3(a) we show a phase diagram of Fo versus 
Fp highlighting the different phases. For weak pinning, 
Fp < 0.06, the depinning transition occurs elastically in 
a single step where all of the vortices begin moving at 
the same time. The effect of removing one pinning row 
is negated by the fact that the vortex ground state for 
weak pinning forms a distorted triangular lattice by shift- 
ing half the vortices off the pinning sites, as shown in the 
inset of Fig. 3(a). The presence of the pin- free channel 
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FIG. 3: The dynamic phase diagram Fo vs Fp for the system 
in Fig. 1. Filled circles: depinning transition. Open circles: 
transition from the SCF to the SB state. Hatched squares: 
transition from the SB state to the moving state. Inset: Vor- 
tex positions (filled circles) and pin positions (open circles) 
in a small portion of the sample illustrating the depinning 
process for Fp < 0.6. Depinning occurs in a single step from 
a partially pinned elastic solid where half of the vortices are 
shifted off the pinning sites, as shown, (b) The phase dia- 
gram Fd vs S for the same system with Fp = 0.14, where 
(5 is a measure of strength dispersion for the pinning sites. 
As 5 increases, the width of the SB region grows. Inset: For 
S — 0.035, there are more steps in {Vx} vs Fd at the onset of 
the shear banding and the large oscillations are replaced with 
a fluctuating or chaotic regime. 



the opening of additional mobile rows of vortices. The 
steps are followed by a transition to a more strongly fluc- 
tuating regime where additional moving rows open and 
close intermittently. 

In summary, we propose a simple periodic pinning ge- 
ometry for vortices in superconductors in which a rich 
variety of complex shear banding phenomena can be re- 
alized. A single row of pinning sites is removed to create 
an easy flow channel. Under an applied drive the vor- 
tices in the channel move first but due to their coupling 
with the vortices in the adjacent pinning sites, vortex 
motion along the neighboring rows is also nucleated. As 
the global depinning transition is approached, bands of 
moving rows of vortices propagate away from and toward 
the channel, creating giant spatio-temporal oscillations 
which can be seen in transport signatures. As disorder 
is added to the periodic pinning array, the velocity-force 
curves develop discrete steps and the oscillating phase is 
replaced by a strongly fluctuating phase. These effects 
are similar to shear banding phenomena in complex fluids 
where different dynamical phases can coexist and couple. 

This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 
DE-AC52-06NA25396. 



merely removes the degeneracy of the ground state and 
ensures that this shift occurs along the cc-direction. For 
Fp > 0.078, the initial depinning occurs plastically and 
the vortex motion first occurs in the pin-free channel. 
The depinning force saturates with increasing Fp since 
the vortices in the pin-free channel are pinned by their 
interactions with the neighboring pinned vortices, and 
this interaction strength does not depend on Fp. The 
transition to the moving phase shifts linearly to higher 
Fd with increasing Fp. Open circles in Fig. 3(a) indicate 
the transition from the SCF to the SB giant oscillation 
regime where multiple rows of vortices are moving. For 
strong pinning, the SB oscillation region is narrow; how- 
ever, the addition of disorder or an increase in the size 
of the pinning sites can extend the width of the SB re- 
gion. To examine the effect of disorder, we changed the 
pinning strength from uniform to a normal distribution 
with mean Fp and width 6. In Fig. 3(b), we plot the 
onset of the different phases for Fjj versus 6 for a system 
with Fp = 0.14. As S increases, the depinning transi- 
tion to the SCF regime is unchanged since the depinning 
force is determined by the fixed vortex- vortex interaction 
strength. The onset of the SB regime shifts to higher Fd 
with increasing disorder. For 6 < 0.035, the dynamics 
are the same as those illustrated in Fig. 1(a); however, 
for larger S the shear banding occurs in a series of steps 
and the large oscillation region is replaced with a strongly 
fluctuating or chaotic regime. In the inset of Fig. 3(b) 
we plot (Vx) vs Fd near the onset of the SB regime for 
6 = 0.035. Here there are four steps that correspond to 
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